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ABSTRACT 

Dynamical  systems  that  are  reversible  in  the  sense  of  Moser  are 
investigated  and  bifurcation  of  trajectories  connecting  saddle  points  from 
stationary  solutions  is  studied.  As  an  application,  reaction-diffusion  models 
in  one  space  dimension  are  considered.  These  equations  are  studied  in  the 
neighborhood  of  a  point,  where  the  set  of  spatially  homogeneous  solutions 
displays  a  Hopf  bifurcation.  It  is  shown  that  from  such  a  point  branches  of 
solutions  bifurcate,  which  can  be  described  as  waves  travelling  to  or  from  a 
center.  These  waves  may  be  exponentially  damped  at  infinity  or  not.  They  can 
be  regarded  as  one-dimensional  analogues  of  "target  patterns"  or  "spiral 
waves" • 
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SIGNIFICANCE  AND  EXPLANATION 

In  experiments  on  reaction-diffusion  systems,  e.g.  the  famous  Belousov- 
Zhabotinskii  reaction,  patterns  of  rotating  spirals  or  propagating  concentric 
rings  are  observed.  These  patterns  have  found  considerable  theoretical 
interest  in  the  recent  literature .  In  this  paper  we  give  a  rigorous  proof  for 
the  existence  of  certain  solutions  to  reaction-diffusion  equations .  The 
qualitative  features  of  these  solutions  are  such  that  they  may  be  regarded  as 
one-dimensional  analogues  of  these  patterns  •  Mathematically,  the  problem  is  a 
degenerate  case  for  bifurcation  leading  to  trajectories  connecting  critical 
points  in  a  reversible  system . 
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The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC ,  and  not  with  the  author  of  this  report. 


BIFURCATION  OF  SINGULAR  SOLUTIONS  IN  REVERSIBLE  SYSTEMS 
AND  APPLICATIONS  TO  REACTION-DIFFUSION  EQUATIONS 

M.  Renardy 

0.  INTRODUCTION 

Recently  Kirchgassner  and  Scheurle  [19],  [20]  have  constructed  bounded 
nonperiodic  solutions  (which  they  call  "singular")  of  reversible  systems  as 
"envelopes"  of  periodic  solutions  with  infinitely  increasing  periods.  Under 
appropriate  hypotheses,  they  prove  that  there  are  branches  of  such  solutions 
bifurcating  from  a  stationary  solution.  In  [33]  I  have  developed  a  different 
approach  to  these  solutions.  In  particular,  it  turned  out  that  the  singular 
solutions  are  in  fact  trajectories  connecting  saddle  points .  The  present 
paper  presents  some  extensions  to  the  results  of  [33]  .  A  global  bifurcation 
theorem  for  singular  solutions  is  shown.  Moreover,  applications  to  reaction- 
diffusion  models  are  discussed.  We  obtain  solutions,  which  may  be  oonsidered 
one-dimensional  analogues  of  the  patterns  of  concentric  rings  or  spiral  waves 
observed  e.g.  in  the  Belousov-zhabotinskii  reaction. 

In  order  to  keep  the  paper  essentially  self-contained,  the  main  results 
of  [33]  are  reviewed  here,  the  reader  will,  however,  be  referred  to  [33]  for 
some  of  the  proofs.  The  paper  is  organized  in  two  parts:  In  the  first  (§§  1- 
5)  we  deal  with  reversible  systems  on  an  "abstract"  level,  in  the  second  part 
(§§  6,7)  we  consider  reaction-diffusion  equations. 

We  study  a  differential  equation  of  the  form 
(0.1)  u  **  —  *  A(U)u  +  B(U,u)  . 

Here  u  is  a  real  parameter,  and  u  lies  in  a  Banach  space  Y.  A(y)  is  a 

linear  operator  in  Y  composed  of  a  bounded  part  depending  smoothly  on  y 

and  an  (in  general  unbounded)  operator  AQ  satisfying  certain  semigroup 

conditions,  and  B  is  a  smooth  bounded  operator  from  R  x  Y  into  Y 

2 

satisfying  #B(y,u)<  *  0(<u<  ).  Our  principal  assumption  is  that  (0.1)  is 
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reversible  in  the  sense  of  Moser  [26] .  This  means  that  there  is  a  linear 

2 

operator  R  e£(Y)  such  that  R  ■  id,  A(U)R  «  -RA(y),  and 
B(M,Ru)  -  -RB(y,u).  Me  assume  that  for  y  <  0  the  spectrum  of  A(y)  has 
positive  distance  from  the  imaginary  axis,  whereas  at  M  *  0  a  pair  of 
eigenvalues  passes  through  0  and  becomes  imaginary  for  y  >  0  .  Under  these 
conditions,  Kirchgassner  and  Scheurle  [19], [20]  have  proved  that,  for  each 
y  in  a  positive  neighborhood  of  0,  there  exists  a  one-parameter  family  of 
periodic  orbits  centered  at  the  origin.  They  also  prove  that,  under  certain 
additional  assumptions,  bounded  nonperiodic  solutions  can  be  constructed  as  a 
limit  of  these  periodic  solutions,  the  convergence  being  uniform  in  bounded 
intervals .  These  limiting  solutions  are  called  "singular"  . 

In  [33]  an  alternative  approach  to  these  solutions  was  given.  The 
existence  proof  is  not  based  on  approximation  by  periodic  orbits t  on  the 
contrary,  it  is  used  that  singular  solutions  are  isolated  in  a  suitable 
function  space.  Parts  of  the  proof  use  ideas  related  to  those  employed  in 
[9],  [21]  and  [28].  A  bifurcation  parameter  e  is  introduced,  leading  to  a 
reduced  equation  for  e  =  0 .  This  reduced  equation  is  not  the  linearization 
of  (0.1),  hut  a  nonlinear  approximation  of  "Ginzburg-Landau"  type  [12].  For 
the  reduced  problem,  singular  solutions  can  be  given  explicitly.  A  refined 
version  of  the  implicit  function  theorem  is  used  to  prove  the  existence  of 
singular  solutions  for  e  0.  Whereas  the  nature  of  the  singular  solutions 
remains  an  open  problem  in  the  work  of  Kirchgassner  and  Scheurle,  it  becomes 
clear  here  that  they  are  in  fact  trajectories  connecting  saddle  points  • 

Under  generic  assumptions  we  have  to  distinguish  two  different  cases, 
which  can  be  shown  to  correspond  to  the  two  different  reversibility  conditions 
in  [19].  In  the  first  case  a  two-sided  branch  of  stationary  solutions  emerges 
from  the  point  u  *  0,  y  ■  0:  For  y  >  0,  the  point  u  *  0,  or  respectively. 


for  M  <  0,  the  bifurcating  fixed  point  are  connected  to  themselves  by  a 
biasymptotic  trajectory.  Using  a  different  method,  this  case  was  also 
investigated  in  [21].  In  the  second  case  a  one-sided  branch  of  stationary 
points  bifurcates  from  0,  i  .e .  either  for  p  >  0  or  for  u  <  0  we  have  two 
new  fixed  points .  it  depends  on  the  direction  of  the  bifurcation  whether 
these  are  saddle  points  or  not.  If  they  are,  they  are  connected  to  each  other 
by  two  trajectories  .  In  both  cases  one  of  the  solutions  represented  by  the 
singular  trajectory  is  symmetric  w.r.  to  R,  i.e.  Ru(-x)  =  u(x).  These 
results  are  reviewed  in  §§  1-4  of  this  paper.  In  §  5  we  use  degree  theory  to 
extend  this  local  bifurcation  theorem  to  a  global  result.  It  is  shown  that 
branches  of  singular  solution  can  terminate  only  by  one  of  the  two  following 
mechanisms : 

1  .  The  asymptotic  stationary  point  loses  the  property  of  being  a 
saddle . 

2.  A  suitably  defined  norm  of  the  solutions  tends  to  infinity. 

There  are  numerous  examples  for  the  occurence  of  singular  solutions  in 

physical  problems.  The  research  of  Kirchgassner  and  Scheurle  was  motivated 
by  the  study  of  stationary  solutions  to  the  Benard  and  Taylor  problems. 

Further  examples  occur  in  classical  mechanics,  in  the  theory  of  water  waves 
[9],  in  the  theory  of  Josephson  junctions  and  in  nonlinear  optics  [3]. 

In  the  second  part  of  this  paper  we  consider  applications  to  reaction- 

diffusion  equations  in  one  space  dimension.  These  equations  have  the  form 

«  • 

(0.2)  Diui  =  f1(y,u1,...,u  )  +  u± 

The  stands  for  the  derivative  with  respect  to  the  space  variable  x,  and 

the  *  stands  for  the  derivative  with  respect  to  time  t.  The  f^  are 

assumed  to  be  smooth  functions.  We  restrict  our  attention  to  solutions 
periodic  in  time.  Putting  u^  =  v^ ,  (0.2)  can  be  rewritten  as  a  system. 


which  is  reversible  under  either  of  the  mappings 

*  fp  m 

R  *  (u^v^)  ♦  *ui'"vi*  or  R  s  <«i(t),vi(t))  ♦  (^(t+j),  -v^t+j)) 

(T  denotes  the  temporal  period).  We  study  solutions  of  (0.2)  in  the 

neighborhood  of  a  point  where  the  system 

(0.3)  f .  (li,u. ,  ...,u  )  +  u  -  0 

ii  n  l 

undergoes  a  Hopf  bifurcation  (this  has  been  shown  to  occur  e.g.  for  the 
Brusselator  [1],  [12]  and  the  Belousov-Zhabotinskil  reaction  [6],  [13],  [17], 
[18],  [27],  [31]).  For  (0.2)  this  leads  to  a  situation,  which  has  analogies 
to  the  one  considered  in  the  first  part  of  the  paper,  but  is  more  complicated) 
the  eigenvalue  0  occuring  in  the  linearized  problem  for  p  »  0  now  has  the 
algebraic  multiplicity  4  and  the  geometric  multiplicity  2  rather  than  2 
and  1  .  Again  we  introduce  a  bifurcation  parameter  £,  and,  taking  into 

account  only  the  terms  of  lowest  order,  we  obtain  a  reduced  equation. 

This  reduced  equation  turns  out  to  be  the  simplest  case  of  a 
"A-co-system" .  A-io-systems  have  been  introduced  as  one  theoretical  concept  to 
explain  patterns  of  concentric  rings  ("target  patterns")  or  rotating  spirals 
("spiral  waves")  occuring  in  chemical  reactions  (see  e.g.  [5],  [8],  [10], 

[11],  [15],  [22],  [23],  [29],  [34],  [35],  [36],  [38]).  Solutions  of 
X-w-systems  were  investigated  in  particular  by  Kopell  and  Howard  [15],  [22], 
[23]  in  one  space  dimension  and  by  Greenberg  [10],  [11]  in  two  space 
dimensions . 

Our  analysis  focusses  on  two  specific  solutions,  which  we  can  give 
explicitly  for  e  =  0.  Again  a  generalized  implicit  function  theorem  is  used 
to  prove  persistence  of  solutions  with  the  same  qualitative  properties  for 
e  7*  0.  The  first  type  of  solutions  are  temporally  periodic,  symmetric  with 
respect  to  R,  and  approach  a  constant  as  x  +  ±»,  asymptotically  they  can  for 
large  |x|  be  described  as  exponentially  damped  waves  propogating  in  opposite 


directions  for  x  positive  and  x  negative •  The  solutions  of  the  seoond 

A 

type  are  temporally  periodic,  symmetric  with  respect  to  R,  in  the  limit 
x  ♦  they  approach  periodic  wave  trains,  and  the  directions  of  propogation 

are  again  opposite. 

The  first  kind  of  solutions  can  be  regarded  as  a  one-dimensional  analogue 
of  target  patterns,  whereas  the  seoond  kind  are  one  dimensional  spiral  waves 
[38]  .  In  the  case  e  ■  0,  the  latter  solutions  coincide  with  one  of  the 
solutions,  for  which  existence  was  proved  fcy  Kopell  and  Howard  in  [23]  .  For 
e  0  these  solutions  were  discussed  on  a  formal  level  by  Cohen, 

Hoppenstaedt  and  Miura  [4]  . 

It  remains  an  open  question  in  general  whether  the  solutions  under  study 
here  can  be  stable  or  not.  In  chapter  7  we  investigate  the  stability  of  the 
solutions  of  the  first  type  for  a  special  range  of  parameter  values  and  find 


they  are  unstable 


I 


BIFURCATION  OP  SINGULAR  SOLUTIONS  IN  REVERSIBLE  SYSTEMS 


1 .  Formulation  of  the  problem 

He  consider  a  differential  equation 

(1.1)  ~  *  u*  ■  A(p)u  +  B(u,u), 

OX 

where  U  is  a  real  parameter  and  u  is  in  a  Banach  space  Y.  We  assume: 

(i)  A (4)  is  of  the  form  A(U)  -  AQ  +  A^  ( M ) ,  where  AQ  A(0)  is  a 

closed,  densely  defined  linear  operator  in  Y  and  A^CtO  @  £(Y)  is  a 
00 

C  -  function  of  U. 

(ii)  B  :  R  x  Y  ♦  Y  is  of  class  C  and  lB(li,u)l  *  0(lul2)  as  u  ♦  0 . 

(iii)  Equation  (1.1)  is  reversible  in  the  sense  of  Moser  [26],  i.e.  there 
exists  a  linear  isometry  R  e  LVt)  such  that 

R2  =»  id,  A(U)R  «  -RA( M)  and  B(P,Ru)  “  -RB(P,u). 

(iv)  Aq  has  an  isolated  algebraically  two-fold  but  geometrically  simple 
eigenvalue  0 . 

Let  N  denote  the  generalized  nullspace  of  Ag  ,  and  M  the 
complementary  subspace  of  Y  which  is  invariant  under  AQ .  It  easily 
follows  from  (iii)  that  M  and  N  are  invariant  under  R.  Moreover, 
it  is  not  difficult  to  prove  that  R|^  has  the  simple  eigenvalues  +1 
and  -1 . 


(V) 

M  has  a  decomposition 

M  =  M+  +  M~,  where 

M+ 

and  M~  are 

invariant  under  A0,  M~» 

■  RM+  Moreover,  -AQ 

1  ♦ 

M 

generates  a  strongly 

continuous  semigroup  of 

negative  type ,  i .e . , 

for 

x  >  0  we  have 

x  — y x 

•  e  o  •  <  Ce  with  positive  constants  C  and  Y  •  It  is  a  simple 

•  +A  x  «Yx 

consequence  that  on  M~  we  have  le  o  I  <  Ce  . 

we  write  u  *  (v,w,z),  where  v  and  w  denote  the  components  in  N  and 

z  e  M.  Without  restricting  generality  we  may  assume  that  R  takes  (v,w) 

to  (v,-w) .  Equation  (1.1)  is  then  rewritten  as  follows: 


3  *  2 

V'  ■  a(y)w  +  Y(y)vw  +  a(y)w  +  wb  (y)z  +  0(|v|  |w|  + 

|w|3|v|  +  »z»(|y|  +  |v|  +  lzl  +  | w | 2 ) ) 

(1.2) 

w'  =  0(y)v  +  <S(M)v2  +  C(y)w2  +  0 ( | v | 3  +  |w|2|v|  +  lzl  • 

(  M  +  |w|  +  9zi  +  | U t > ), 

z*  =  A (U)z  +  w2a(y)  +  0(*zi(|v(  +  | w |  +  9zl)  +|v|2  +  |v|  |w| 

+  | w | 3  +  |y | ( |v|  +  |w| ) ) , 

where  a(y),  0(y),  5(y),  ?(y),  and  a(y)  are  real  numbers, 

*  *  -v 

a(y)  e  M,  b  (u)  e  M  (the  dual  of  M) .  and  A(y)  is  a  linear  operator  in 
M.  We  shall  distinguish  the  following  generic  cases: 

Case  1 : 

0Q  =  0(0)  =  0,  01=  0(y)ly=o  ?  o,  sq  =  6(0)  ft  o. 

2 

Then  we  put  U  =  ±g  and  introduce  the  scaling 
v  -*■  e2v,  w  *  e3w,  z  +  g3z,  x  ♦  ^  .  We  obtain: 


V 

=  aQw  +  0(| g | ) 

(1.3) 

« 

w 

=  ±0^  +  6qV2  +  0  ( |  e  | ) 

« 

ez 

=  A(0)z  +  0 ( | e I ) 

Case  2 : 

a  =  a(0)  =  0,  a  =  a(y) |  *  0,  0(0)  *  0,  0(O)a(O)  -  Y(O)C(0)  -  0, 

o  i  ay  y=o 

where  a{0)  =  a(0)  +  b*(0 ) A(0 ) -1a(0) . 

^  2'-  - 1 

We  then  put  z  =  z-w  A(0)  a(0)  and  introduce  the  scaling 

2  %  2 N  x 

v  ♦  g  v,  w  +  Gw,  z+Gz,x+— .  We  obtain 

v  =  +  Yqvw  +  <7gW  +  0(*z*  +  |g|), 

(1.4)  w  =  0Qv  +  CQw2  +  0 ( | G | ) , 

X  •  V  X 

Gz  =  A(0)z  +  0(  |e()  , 


-7- 


I 
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2.  Singular  solutions  for  e  =  0 


We  are  now  going  to  discuss  the  existence  of  trajectories  connecting 
saddle  points  for  e  =  0.  Since  A(0)  is  nonsingular ,  e  »  0  immediately 
yields  z  =  0  both  for  (1.3)  and  (1.4),  and  we  are  left  with  a  two  dimensional 
problem  in  either  case.  We  start  with  the  easier  case  1. 

Case  1 : 

We  assume  >  0,  8^  <  0,  which  can  easily  be  achieved  by  replacing 
Vi  and  v  by  -U  and/or  -v  if  necessary.  For  e  *  0  (1.3)  reads 


(2.1) 


v  =  aQw, 

i  2 

w  =  v  , 


which  is  a  Hamiltonian  system,  i.e.  the  Hamiltonian 


8  26 

„  2  -  1  2  ^  0  3 

H  =  W  +  —  V  -  T—  V 

%  3o0 


is  constant  along  trajectories.  This  implies  the  following 
Proposition  2.1: 

If  the  plus  sign  is  chosen  in  (2.1)  (corresponding  to  y  >  0),  then  the 
fixed  point  v  “  -8^/6^,  w  =  0  is  a  saddle  point,  which  is  connected  to 
itself  by  a  separatix,  if  the  minus  sign  is  chosen,  the  same  holds  for  the 
fixed  point  0 . 


Case  2: 


(1.4)  reads  for  e  «  Os 


(2.2) 


3 

±<*|W  +  1Qvw  +  oQw  , 

v +  V2  • 


Again  we  may  assume  8Q  >  0,  <  0.  Stationary  solutions  of  (2.2)  are  given 


=  w  =  0  or  v  =  ±Vo/(BoV  Vo,  ,  w  =  ;  a^Vo  -  W 


Hence  nontrivial  fixed  points  exist  for  U  >  0  (i.e.  for  the  choice  of  the 

4-  sign  in  (2.2))  if 


(2.3) 


Vo  "  Vo 


>  0. 


A  simple  calculation  shows  that  in  this  case  the  nontrivial  fixed  points 
are  saddle  points.  Separatrices  connecting  the  two  nontrivial  fixed  points 
are  found  as  follows: 

2 

We  look  for  invariant  parabolae  of  the  form  v  =  aw  +  b.  This  curve  is 
invariant  under  the  flow  of  the  differential  equation  (2.2)  iff: 

I  f 

v  =  2aww  .  If  (2.2)  is  inserted  into  this  equation,  a  short  calculation 
shows  that  there  are  in  fact  two  invariant  parabolae  for  (2.2),  namely,  we  get 


a 


1  ^  /  l  2  ^ 

^  Y  -  ±  /(-  yn-  +  2a  0 

2  0  0  »  2  0  0 _ 0  C 

20. 


t 


±a. 


b  = 


2  Y0  “  ?0  * 


i/'F'iT 


v 


+  2<,o  Bo 


It  can  be  checked  that  the  nontrivial  fixed  points  are  in  fact  on  these 
parabolae.  We  leave  the  calculations  to  the  reader.  That  the  two  parabolae 
are  the  only  trajectories  connecting  the  two  saddle  points  follows  from  the 
uniqueness  of  stable  and  unstable  manifolds,  for  which  we  refer  to  [30]  . 
Altogether  we  find 
Proposition  2,2: 

If  (2.3)  is  satisfied  and  the  +  sign  is  chosen  in  (2.2)  (which 
corresponds  to  P  >  0  ),  then  there  exist  two  saddle  points,  whidi  are 


connected  to  each  other  by  two  trajectories.  The  solutions  of  (2.1)  and  (2.2) 


3 .  The  linearization  at  the  singular  solutions 


All  the  trajectories  aonnecting  saddle  points  which  we  have  found  in  §  2 
are  symmetric  with  respect  to  the  w-axis .  This  means  that  among  the  one- 
parameter  family  of  solutions  represented  by  such  a  trajectory  there  is  one 
solution  yQ(x)  =  (vQ(x),  wQ(x))  satisfying  RyQ(x)  *  yQ(-x).  Me  denote  the 
linearization  of  (2.1)  resp.  (2.2)  at  y0(x)  by 

I 

y  *  CQ(x)y 

Using  the  reversibility  and  the  fact  that  RyQ(x)  **  yQ(-x)  one  finds: 

Cq(-x)R  =  -RCq(x).  We  shall  prove 
Theorem  3.1 

For  each  f(x)  *  (f  (x) ,f „(x) )  e  F  *  (f  6  C®(R,R2) |Rf (x)  »  -f(x), 

12  m  d 

/i.\  /  v  \ 

lyj  f'  ' (x) ,  l^m^  *  7 (*)  exist  for  0<k<m)  there  exists  one  and  only  one 
y(x)  =  (v(x),  w ( x ) )  e  Um+1  «  (y  e  c“+1 (R,R2) |Ry(x)  =  y(-x),  Ijl  V<k) (x) . 
y  (x)  exist  for  0<k<m+1 }  solving  the  inhomogeneous  equation 

(3.1)  y’  -  CQ(x)y  =  f 

in  2  2 

Here  C.(R,R  )  denotes  the  Banach  space  of  all  functions  R  ♦  R  having  m 

D 

continuous  bounded  derivatives . 

Proof : 

For  x  +  ±®»,  yQ(x)  converges  to  a  saddle  point,  and  one  easily  concludes 
from  the  stable  manifold  theorem  [30]  that  the  convergence  is  exponential. 
Repeated  differentiation  of  (2.1)  resp.  (2.2)  then  yields  the  result  that  all 
derivatives  of  y^  converge  to  0  exponentially.  This  implies  the  following 
properties  of  CQ(x):  For  x  ♦  ±“,  CQ(x)  converges  to  the  linearization  at  a 
saddle  point,  i.e.  l^ij  Cg(x)  exists,  and  this  matrix  has  one  positive  and 
one  negative  eigenvalue.  Moreover  dm/dxm (C^ (x) )  converges  to  0 
exponentially  for  each  m  >  0.  From  this  we  see  that  it  is  sufficient  to  prove 
the  theorem  for  m  =  0,  the  rest  following  from  (3.1)  by  repeated 


I 


1 


differentiation.  We  rewrite  (3.1)  in  the  form 

I  A  % 

y  -  Cy  -  CQ(x)y  =  f  , 

A  ^  A 

where  c  is  constant,  and  lim  n  (x)  =0.  From  the  fact  that  C  has  one 

x**  0 

positive  and  one  negative  eigenvalue  we  conclude  that  for  each 


f  e  C,  ([X,00),!*  )  the  equation 
b 

(3.2) 


y  -  Cy  =  f 


1  2 

has  a  solution  y  e  C,  ([X,«),R  ),  which  is  determined  up  to  one  arbitrary 

D 

initial  condition  at  x  =  X.  If  1J$  f(x)  exists,  then  l^jj  y(x)  exists  as 

well.  Let  P  denote  a  projection  of  R  onto  a  one-dimensional  subspace  such 

that  the  prescription  of  Py(X)  determines  the  solution  of  (3.2)  uniquely. 

Then  the  mapping  y  ♦  (y-  Cy,  Py(X))  is  an  isomorphism  from 

Ux  =  (y  e  Cfa([X,®),R  )  1 1^53  y(x),  l^gj  y  (x)  exist)  onto 

(F  -  {f  e  CJ[X,”),R2)|  lim  f(x)  exists  })  x  R.  Since  (3.2)  is 
X  D  X** 

autonomous,  the  norm  of  this  isomorphism  and  its  inverse  are  independent  of 
X.  On  the  other  hand,  the  norm  of  the  mapping  y  ♦  CQ(x)y  from  Ux  onto  Fx 
tends  to  0  as  X  ♦  ®,  i .e .  for  sufficiently  large  X  the  mapping 

•  A  x 

y  ♦  (y  -  CyQ-  C  (x)y,  Py(X))  is  still  an  isomorphism  from  Ux  onto  Fx  x  R. 
That  means,  given  any  f  e  FQ,  there  exists  a  bounded  solution  y  to  (3.1) 
on  the  interval  (X,®),  which  is  determined  up  to  one  initial  condition  at  x 
=  X.  Since  the  initial  value  problem  for  (3.1)  is  uniquely  solvable,  this 
solution  on  (X,®)  extends  to  a  solution  on  [0,®),  and  we  have  one  free 
initial  condition  at  x  =  0.  This  remaining  free  initial  condition  is  matched 
by  adding  a  multiple  of  the  solution  yg(x)  of  homogeneous  problem.  The 
requirement  y(x)  ■  Ry(-x)  implies  w(0)  *  0,  which  is  achieved  by  one  and 

I 

only  one  choice  of  the  initial  condition,  since  ( 0 )  /  0.  This  determines 
the  "half-sided"  solution  on  (0,®)  uniquely.  Since  Rf(x)  =  -f(-x),  we 
have  further:  If  y(x)  solves  (3.1)  on  [0,®),  then  Ry(-x)  solves  (3.1)  on 


1 


(-•,0]  (use  the  fact  that  CQ  (-x)R— Rcq  (x) ) .  Hence  the  half-sided  solution 
extends  to  a  solution  on  all  of  R,  and  the  theorem  is  proved • 


I 


4 .  Existence  of  singular  solutions  for  e  ^  0 


The  goal  of  this  paragraph  is  to  prove  that  there  exists  a  branch  of 
singular  solutions  in  a  neighborhood  of  e  *  0  .  For  this  let  y^(x)  be  as 
in  §  3  and  put  h(x)  =  (v(x),w(x))  -  yQ(x)  .  (1.3)  or  (1.4)  then  takes  the  form 

h  =  CQ(x)h  +  f(e,h,z,x)  , 

(4.1) 

•  % 

£z  =  A(0)z  +  g(e ,h,z,x)  • 

2 

Here  CQ(x)  is  as  in  §  3  and  we  have  lf(e,h,z,x)l  *  0(|e|  +  izl  +  Ihl  ), 

<g(e,h,z,x)<  =  0(|e|).  For  each  m  e  N  the  mapping  (e,h,z)  ♦  (f(e,h,z,x), 

00 

g(e,h,z,x))  is  a  C  -mapping  from  R  x  U  (Y)  to  F  (Y),  where  we  denote 

m  m 

Um<Y)  =  {u  e  C^(R,Y)  |l^jg  u{k,(x)  ,  l^m^  u(k,(x)  exist  for  0  <  k  <  m, 

Ru(x)  =  u(-x) }  and 

Fm(Y)  =  {f  e  C™(R.Y) iLJg  f (k) (x) ,  l^m.  f(k)(x)  exist  for  0  <  k  <  m, 

Rf (x)  -  — f ( x) } 


we  write  (4.1)  as  follows: 

h  -  (^  -  Cgfx))”1  f  (e,h,z,x)  =*  0 


(4.2) 


z  -  (e  ^  -  A(0 ) )  1  g(e,h,z,x)  =  0 

we  wish  to  apply  an  implicit  function  argument  to  establish  the  existence  of 

d 

solutions  in  Um(Y)  for  e  ^  0.  A  problem  is  caused  by  the  fact  that  e  — 

■>. 

is  a  relatively  unbounded  perturbation  to  A (0).  It  is,  however,  not 
difficult  to  show  the  following  (the  details  are  in  [33]): 

Lemma  4.1: 

d  »  -1 

(i)  The  mapping  (£,z)  (£  ~  -  A(0))  z  from  R  x  F  (M)  into 

ax  m 

U  (M)  is  continuous  near  e  =  0  (M  as  defined  in  (iv)  of  §  1) 
in 

d  ^  -1 

(ii)  The  operator  norm  of  ( £  - -  A(0))  :  F  (M)  ♦  U  (M)  is 

ax  m  m 

uniformly  bounded  in  a  neighborhood  of  e  ■  0  . 
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(iii)  The  mapping  (e,s)  ♦  (e  4~  -  AfO))”^  from  R  x  F  (M)  into 

ax  m 

U  .  (M)  it  of  class  Ck. 
m-K 

These  properties  permit  the  use  of  the  following  abstract  theorems  which 
we  proved  in  [32] ,  [33]  . 

Theorem  4.2 

Let  X,Y  and  z  be  Banach  spaces,  u  a  neighborhood  of  (0,0)  in 
X  x  Y,  and  Ft  u  +  z  a  mapping  having  the  following  properties: 

(i)  F(0,0)  -  0. 

(ii)  F  is  continuous. 

(iii)  F  is  continuously  differentiable  with  respect  to  y  for  each 
fixed  x. 

(iv)  D^F(0,0):  Y  ♦  Z  is  an  isomorphism. 

(v)  DyF  is  continuous  at  the  point  (0,0). 

Then  the  equation  F(x,y)  *  0  has  a  unique  continuous  resolution  y  -  f(x) 
in  some  neighborhood  of  (0,0). 

Theorem  4.3 

(k)  (k) 

Let  Y  resp.  Z  (k  »  0,1,...  N)  be  two  hierarchies  of  Banach  spaces 
(k)  (k+1 )  (k)  (k+1 ) 

such  that  Y  Cy  ,  z  C  z  ,  the  imbeddings  being  continuous.  Let 
X  be  a  finite  dimensional  Banach  space  and  F  a  mapping  from  a 
neighborhood  U  of  0  in  X  x  Y^  into  Z^  having  the  following 
properties : 

(i)  F(U  O  (X  x  Y(k)))  C  z(k)  k  -  0,1,...  N 

(ii)  For  each  fixed  k,  F^s*  F|y  n  (x  x  Y^)  satisfies  the  conditions 

(k) 

of  theorem  4.2,  when  it  is  considered  as  a  mapping  from  X  x  Y  into 
(k) 

Z  .  For  x  fixed,  F^(x, .)  is  a  smooth  (i.e.  sufficiently  often 
differentiable)  mapping. 

(iii)  F:  X  x  Y(k)  ♦  Z(k+w)  is  of  class  Cm  for  each  k  -  0,  1,...,  N 
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and  m  <  N  -  k 


(iv)  The  mapping  (x^u* , . .  .,u^  )  +  z  -  D  i  j  F(x,y)  (u1 , . . . ,u"* )  is  a 

x  y 

continuous  mapping  from  X  x  Y^  x  (Y^)^  into  £*(X,Y  )  . 

Then  the  following  holds: 

The  solution  y  *  f(x)  e  Y^  existing  by  theorem  4.2  is  a  (^“-function  of 
x  in  some  neighborhood  V  of  0,  if  y  is  regarded  as  an  element  of 

Rl 

»<mi. 

(V)  (\c) 

Identifying  X  with  R,  Y  and  z  with  u  .  (Y),  we  get  from 

m-K 

these  theorems . 

Theorem  4.4: 

For  each  m  e  N  there  exists  a  neighborhood  of  (0,0)  e  R  x  U  (Y) 

m 

such  that  in  equation  (4.2)  has  a  unique  resolution 

h  =  h(e),  z  =  z(e)  .  If  this  solution  is  considered  lying  in  u  (Y),  then 

m-k 

u 

in  some  neighborhood  of  e  ■  0,  it  is  a  C  -function  of  e  . 

Remarks 

1 .  In  order  to  carry  out  the  iteration  procedure  in  the  proof  of  the  implicit 
function  theorem,  equations  of  type  (3.1)  have  to  be  solved.  Although  the 
explicit  solution  of  (3.1)  has  not  been  employed  in  the  proofs,  it  can 
easily  be  obtained  modulo  integrations,  since  one  integral  of  the 
homogeneous  equation  is  known . 

2.  in  the  case  1  of  §  1,2  there  are  singular  solutions  approaching  0  for 

x  ♦  ±*  .  Using  the  stable  manifold  theorem  it  can  be  shown  that  the 

convergence  is  exponential,  and  the  singular  solutions  are  therefore  in 

Lp  for  each  p  >  1  .  But  do  they  bifurcate  in  LP?  To  answer  this 

question,  we  must  see  how  the  scaling  introduced  in  §  1  affects  the  LP- 

2 

norm.  In  the  generic  case  considered  here  u  has  been  scaled  by  e 
and  x  by  e  which  gives  a  factor  of  e2  1//p  in  the  LP-norm. 
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Therefore  we  have  a  bifurcation  in  the  space  iP  for  each  p  >  1 .  This 
need  no  longer  be  true  if  degeneracies  occur  and  different  scaling  factors 
roust  be  used.  KQpper  and  Riemer  [24]  have  considered  the  example 

-u  -  |u|rf(u)  *  Xu 

where  f(0)  <  0  and  r  <  1. 

Our  method  described  above  applies  to  this  example  with  the  scaling 
X  =  -e  ,  u  ♦  e  7  u,  x  ♦  e  x  «  In  the  LP-norm  this  gives  a  factor  of 
£2/(r  1)  1/p  '  This  exponent  is  greater  than  zero  if  r  <  2p  +  1 .  This 
agrees  with  the  result  obtained  in  [24]  for  the  case  p  =  2. 
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5 .  Global  existence  of  singular  solutions 


4 


In  the  preceding  chapters  we  have  proved  the  existence  of  branches  of 

singular  solutions  in  the  neighborhood  of  some  bifurcation  point.  This 

chapter  deals  with  the  problem,  how  far  these  branches  can  be  continued.  The 

main  tool  of  the  analysis  will  be  the  theory  of  degrees  of  mappings.  Under 

appropriate  conditions,  a  degree  can  be  associated  with  our  singular  solutions 

in  a  quite  similar  way  as  with  solutions  representing  travelling  waves  [39]  . 

Since  the  definition  of  a  degree  requires  a  compactness  assumption,  we 

impose  the  following  condition  in  addition  to  (i)  -  (v)  of  §  1 s 

(vi)  -A  |  and  A  |  ,  respectively,  generate  analytic  semigroups  and 

0  M  0  m“ 

have  compact  resolvents. 

a  a 

As  a  consequence,  fractional  powers  (-A^)  And  are  defined  as 

operators  in  M~  and  M+,  respectively.  We  say  briefly  that  y  e  ¥  is  on 

a  a  +  a 

D(aq),  if  its  M  -component  is  in  D ( -AQ )  and  its  M  -component  is  in  D(AQ)  . 

It  follows  from  the  compactness  of  the  resolvent  that,  for  any 
a 

0  <  a  <  1,  d(Aq)  is  compactly  embedded  in  Y. 

The  result  we  want  to  show  is  that  branches  of  solutions  can  be  continued 
unless  some  norm  of  the  solution  approaches  infinity.  Obviously,  this  can 
only  be  expected,  if  the  nonlinear  terms  in  the  equation  remain  bounded  on 
bounded  sets .  We  therefore  assume 

(vii)  The  nonlinear  operator  B  maps  bounded  sets  into  bounded  sets. 

In  our  exposition  we  focus  on  the  case  1  of  §1  .  Again  we  assume 
°0  >  <  O’  In  thi®  esse  we  have  proved  that  a  branch  of  trajectories 

connecting  the  saddle  point  0  to  itself  exists  for  V  <  0. 

We  shall  be  concerned  with  solutions  lying  in  the  following  spaces. 
Definition  5.1; 

Let  Y  be  the  same  Banach  space  as  in  §  1 .  Then 
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x”(Y)  {y(x)  e  Cn(R,Y)  jRy(x)  -  y(-x) ,  sup  e° ^ ly ^ (x) i  <  •  for  k  -  0,1,... n) 

(k) 

y  denotes  the  kth  derivative  of  y. 

In  §§  1-4  we  have  shown  that,  for  e  small,  (1.3)  has  two  solutions  in 
x“  (Y)  ,  namely  0  and  the  singular  solution.  Both  are  isolated,  and  the 
linearization  of  (1.3)  at  either  solution  is  nondegenerate. 

For  convenience,  we  let  y  denote  (v,w,z),  L(e)  the  linearization  of 
the  right  hand  side  at  0,  N(e,y)  the  nonlinear  part  of  the  right  hand  side 
of  (1.3),  and  <pE  the  mapping  (v,w,£z)  .  The  (1.3)  can  be  rewritten  in  the 
form 

(5.1)  y  *  { dx  " 

Let  T  m  f c ^  ,c^]  denote  an  interval  with  the  following  properties: 

1.  Z2>  €1>  Gi  s®**! 

2.  For  each  e  e  T,  the  solution  0  of  equation  (1.3)  is  a  saddle  point, 
and  the  real  parts  of  the  eigenvalues  of  'l(s)  have  a  positive 
distance  5  (uniformly  in  e  )  from  the  strip 

(A  e  C|  -o  <  Re  X  <  o} 

The  crucial  property  for  a  global  bifurcation  result  is  the  following 
Propositon  5.1 

Let  D  be  any  ball  in  X^(Y)  .  Then  the  right  hand  side  of  equation 
(5.1)  represents  a  completely  continuous  mapping  from  T  x  D  into  x”(Y) 

Proof : 

It  follows  from  assumption  (vii)  that  N  is  continuous  and  bounded  into 

X2o  ^ *  *  Moreover  *  our  assumptions  on  and  r  imply  that 

d  -1  -1 

(—  -w  L(e))  is  a  bounded  continuous  mapping  from 
ax  £ 

r  x  X^fY)  into  XP(D(A>  for  some  a,n  >  0  and  C  <  min  (2 c,a+i).  The 
rest  follows  from  the  compactness  of  the  embedding  D(A^)  *  Y  and  the  Arzela- 
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Asooli  theorem 


Degree  theory  now  gives  the  following  result  (cf.  [37]): 

Theorem  5.2: 

The  branch  of  singular  solutions  to  (1.3)  provided  by  theorem  4.4  must 
either  leave  T  or  the  norms  of  the  solutions  are  unbounded  in  x”(Y). 
Remarks : 

1.  The  reader  is  cautioned  that  unboundedness  in  (Y )  does  not  mean 

unboundedness  in  the  sup-norm.  In  fact,  branches  of  singular  solutions 
can  terminate  by  a  mechanism  like  the  one  indicated  in  the  following 
sequence  of  diagrams . 


For  any  o  >  0,  the  norm  in  x”  tends  to  09  when  passing  from  situation 
(a)  to  situation  (b),  and  the  branch  of  solutions  connecting  0  to  itself 
terminates  its  existence.  Similar  phenomena  must  be  expected,  when  singular 
solutions  come  close  to  an  invariant  set  other  than  a  stationary  point. 


2.  The  above  theorem  does  not  exclude  the  possibility  that  a  branch  goes  back 
to  e  «  0.  However,  branches  can  not  return  to  the  point  they  bifurcated 
from,  as  one  can  see  from  a  local  uniqueness  result.  Namely,  there  exists 
a  center  manifold  z  =»  g(v,w,u)  for  equation  (1.1),  and  for  w  near  0 
all  uniformly  small  solutions  must  lie  on  the  center  manifold.  This 
follows  from  the  following  argument: 

A  A 

If  we  put  z  *  z  -  g(v,w,W),  we  obtain  for  z  a  differential 
equation  having  the  following  form 

*  I  ^  A  A  A  * 

z  -  A(0)z  +  0( ( | u |  +  | v |  +  |w|))  Izl)  +  0(lzl  ) 

A 

For  small  |v|,  |w| ,  |w|  ,  the  only  small  solution  of  this  is  z  =  0,  as 
follows  from  the  implicit  function  theorem.  From  this  and  elementary 
considerations  about  two-dimensional  flows  it  follows  that  near  V*  =  0 
the  saddle  point  0  and  the  singular  solution  provided  by  theorem  4.4  are 
the  only  small  solutions  in  . 

Ftor  the  case  2  of  §  1,  analogous  considerations  are  possible,  but  now 
L(e)  should  be  replaced  by  an  appropriate  operator  L(e,x)  which 
converges  to  the  linearization  at  the  limiting  fixed  points  for  x  ♦  ±®. 
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II.  APPLICATIONS  TO  REACTION-DIFFUSION  MODELS 


6 .  Oscillating  sincrvlar  solutions  connected  with  Hopf  bifurcations 


We  consider  a  general  chemical  reaction  model  given  by  an  equation 

3u 


32u 


(6.1 )  ^  =  F( M,u)  +  D 

3x 

where  u  e  R  ,  M  e  R.  F  is  a  smooth  nonlinear  function.  D  is  a  diagonal 
matrix,  which  is  strictly  positive  definite.  We  assume 

(i)  For  M  in  some  neighborhood  of  0,  there  exists  a 

solution  u  =  Ug(M)  e  Rn  to  the  equation  F(U,u)  =  0  ,  and  uQ 

CO 

is  a  C  -function  of  M 

(ii)  The  matrix  D^FfOjUgfO))  has  the  algebraically  simple  imaginary 
eigenvalues  ±iwQ ,  the  rest  of  the  spectrum  lies  in  the  left  half 
plane . 


(iii)  Let  1(U)  denote  the  branch  of  eigenvalues  of  D^F( M,uQ (M) ) , 

which  goes  through  iw_  at  V*  =  0  .  Then  ■§-  Re  X(u)|  *  0. 

0  dM  M=0 

(iv)  For  each  Y  >  0,  the  spectrum  of  D^FfOjUgfO))  -  DY  lies  in  the 
left  half  plane  . 

It  is  well  known  [14]  that  conditions  (i)-(iii)  imply  the  existence  of  a 
branch  of  x-independent  time-periodic  solutions  emerging  from  the  bifurcation 
point  u  =  uQ(0),  M  =  0.  The  conditions  (i)-(iv)  have  been  verified  in  quite 
a  few  reaction  diffusion  models,  e.g.  the  "Brusselator"  [1],  [12]  and  the 
Field-Noyes  model  of  the  B-Z  reaction  [6],  [13],  [17],  [18],  [27],  [31]. 

In  this  chapter  we  shall  study  time-periodic  space- depen dent  solutions  of 
equations  (6.1).  We  rewrite  this  equation  in  the  form 


(6.2) 


1  "  -1 c 3u  1  ,  ^ 
-  u  =  D  -  -  F(U,u)J 


The  factor  o>  >  0  has  been  introduced  in  order  to  normalize  the  period  to 
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I 


2*.  After  appropriate  scaling  of  x  we  may  drop  —  on  the  left  side. 
Moreover,  we  write  (6.2)  as  a  first  order  system. 


(6.3) 


•  -1 /3u  1  .  -i 

V  “  D  lit  "  M  F(W'U)J 


Let  now  fcJ(Rn)  denote  the  space  of  all  2i»-per iodic  functions  ys  R  ♦  Rn, 


r  Cl  (]r)  fie) 

for  which  *y*a  :*  L  (|k|  +  1)  ly  I  <  ",  where  the  y  are  the  Fourier 

kSS  ^j^t 

coefficients:  y(t)  »  I  y  e  .We  shall  seek  solutions  to  (6.3)  in  the 

kes 

space  Y  -  {(u,v)|u  e  (Rn)  ,  v  e  (Rn)},  where  m  is  an  arbitrary 

m  m  m— 

positive  integer.  Clearly,  the  mapping  (M,w, (u,v))  (0,-D  ^(u.u))  is 

smooth  from  R  x  (R  \  (o) )  x  Y  into  Y  ,  moreover,  the  operator 

m  m 

9U 

(u,v)  (v,D  -g^)  is  densely  defined  in  Y^  and  closed.  Equation  (6.3)  is 
reversible  under  both  mappings  R  :  (u,v)  **■  (u,-v)  and 

A 

R  :  (u(t)  ,v(t) )  **•  (u(t+ir),-v(t+ir)). 

We  now  discuss  the  spectrum  of  the  linearization  of  the  right  side  of 
(6.3)  at  the  point  to  *  <dQ  ,  y  =  0,  u  *  uQ  (0)  .  The  relevant  properties  of 
this  spectrum  are  described  by  the  following  lemma. 

Lemma  6.1; 

Assume  (i)-(iv)  and  (v)  stated  below  hold.  Then  the  operator 

A  :  (u,v)  **■  (v,D  T-  -  (Dttf  )  D  F(0,uft(U))u)  has  the  Isolated 
« t  0  u  0 

algebraically  four-fold  and  geometrically  two-fold  eigenvalue  0.  Let  N 
denote  the  generalized  nullspace  and  M  a  complementary  invariant  subspace. 
Then  A|„  satisfies  oondition  (v)  of  §  1  (even  (vi)  of  §  5,  but  we  shall  not 
use  this.) 


* 


Proof: 


We  first  note  that  A  acts  Fourier-coraponentwise,  i.e.  we  have 
A(I(u(k),v(1°)eikt  -  I  A(k)'-(k)  -(k)‘-lkt 


where 


(-’  -  1 

\D  ik  -  (DU0)  DuF(0,uq(0))  0 

(k) 

Thus  the  eigenvalues  of  A  are  the  square  roots  of  the  eigenvalues  of  the 
matrix 

D_1ik  -  (D<o0)’1DuF(O,u0(O)) 

(k) 

and  A  has  an  imaginary  eigenvalue  iff  this  matrix  has  a  negative  real 
eigenvalue .  Let  now  be 

(D-1  ik  -  (D<oo)“1DuF(0#u0(0)))y  «  -Yy  Y  >  0 

This  yields 

DuF(°'u0<0))y  “  V1*  “  i5tWoy 

For  Y  *  0  this  is  impossible  by  condition  (iv),  and  Y  ■  0  yields  k  =  ±1 , 

since  ±i<^Q  are  the  only  imaginary  eigenvalues  of  d^f(0 ,uQ (0 ) ) .  we  see 
(k) 

therefore  that  A  has  no  imaginary  eigenvalues  for  k  *  ±1 ,  and  the  only 
imaginary  eigenvalue  for  k  =*  ±1  is  equal  to  0 .  To  prove  the  statement 
concerning  the  multiplicity,  it  must  be  shown  that  X  =  0  is  an  algebraically 
simple  eigenvalue  of 

D_1i  -  (DM^D  F<0,u  (0>) 

It  is  easy  to  see  that  the  nullspace  of  this  matrix  is  spanned  by  the 
eigenvector  of  DuF(0,uQ(0))  to  the  eigenvalue  ,  whence  the  eigenvalue 
X  *  0  is  geometrically  simple.  Assume  now  that 
DuF(0,uQ(0))y  -  iwQy 

and 

y  *  (D  \  -  (DW^)  1DuF(0,Uq(0 ) )z 

This  yields 

Dy  >=  (iu>  -  D  F(0,u  (0)))z 
0  u  u  o 
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We  shall  therefore  assume 

(v)  With  y  denoting  the  eigenvector  of  DuF(0,Ug(0))  corresponding  to  the 
eigenvalue  iw^,  °y  is  not  in  the  range  of  iwQ  -  DuF(0,u0(0)) . 

Clearly  (v)  is  a  generic  condition,  which  guarantees  that  the  eigenvalue  is 
algebraically  simple  (in  particular,  (v)  follows  from  (ii)  if  D  is  the 
identity  matrix)  . 

We  must  now  verify  that  the  spectrum  of  A  is  actually  given  by  the 
(k) 

eigenvalues  of  the  A  and  that  (v)  of  §  1  holds.  This  will  be  a 
consequence  of  the  followings 

There  exists  an  isomorphism  T  of  acting  Fourier  componentwise: 

Y  (k)  (k),  ikt  r  _(k)  (k)  (k).  lkt 

T  I  (u  ,v  )e  -  l  T  (u  ,v  )e 

kez  kez 

such  that  for  large  |k|,  let  us  say  for  |k|  >  kQ  the  matrix 
(T*k*)  ^A^T^  consists  of  a  diagonal  part  and  a  rest  term  which  has  a 

-  Vi 

norm  of  the  order  of  magnitude  |k|  *  .  Namely,  put  for  |k|  >  kQ 
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Odkl"1) 


(T  )  A  T 


-VjF4  0(|ltl*,/2) 


0(1) 


\ 


lf£.  I 


(k)  -1  (k)  (k) 

i .e .  (T  )  A  T  is  the  sum  of  diagonal  matrix  and  a  remainder  term 

_  1/_ 

which  has  an  operator  norm  of  the  order  of  magnitude  |k|  2  (recall  the 

definition  of  the  norm  in  Y  ) . 

m 

It  is  now  a  simple  consequence  of  the  Hille-Yosida  theorem  [25]  that 
property  (v)  of  §  1  holds  for  the  diagonal  part,  since  the  necessary  resolvent 
estimates  are  in  this  case  trivial,  A  perturbation  argument  shows  that  the 
same  holds  for  the  full  operator  T  ^AT,  and  hence  for  A,  This  concludes 
the  proof  of  the  lemma . 

Again  y  shall  denote  the  eigenvector  of  DuF(0,uQ(0))  corresponding  to 
the  eigenvalue  iU)Q  •  We  decompose  (u,v)  as  follows 


(u , v )  =  (u0(ij),0)  +  ai(y,0)eit  +  a^{y,0)e  +  a2(0,y)eit;  +  <*2  (0,y)e  +  z 


yhere  '0t2  e  C  ar*d  z  e  M  (as  defined  under  lemma  6.1),  Equations  (6,3) 
then  assume  the  form 


°1 

aj  -  “Vi  4  "2V1  4  V'V  4  C2U'V  4  <“■'  -  %’>Vi  4  •• 

MO),  +  a*a,  *  t^a3  1  ... 

Here  a^,  a2»a3  are  complex  numbers,  C1  2  :  M  x  C  +  C  are  bilinear 

■V 

operators  and  the  d  are  vectors  in  M,  A(0)  denotes  A|  .  The  dots 

1  M 


indicate  higher  order  terms . 


Analogously  as  in  §  1  we  put  z  =  z  -  A(0)  {o^d^  +  °.|ajd2  +  o^d^}  and 

introduce  the  scaling 

0^  eai ,  a2  £2a2,  M  =  ±e2,z  ♦  e2z,  <»>  ♦  e^u>,  x  ♦  e  ^x.  We  then 


-2 
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obtain  an  equation  of  the  following  form 


(6.4) 


aj  -  ±  aial  +  a2a^o1  +  a3u«1  +  0(|e|  +  «zl) 
ex'  -  A(0 )z  +  0( |e|) 


For  e  »  0  this  reduces  to 


(6.5) 


“l  “ 


+  a  <* 
*  1  1 


a3°*1 


'  2— 
a2°tai 


Apart  from  a  factor  of  v  -  (where  b*  denotes  the  left  handed 

(b* ,Dy) 

eigenvector  of  DuF(0,uQ(0))  corresponding  to  iu>Q,  normalized  such  that 

(b*,y)  =1),  the  coefficients  a^,  a2  and  a^  are  the  same  which  determine 

3 

the  Hopf  bifurcation  at  the  order  e  .  it  follows  from  assumption  (iii)  that 

Re  —  is  negative.  Moreover,  a3  is  equal  to  iv.  id  is  an  unknown 

variable,  which  has  to  be  determined. 

Equation  (6.5)  is  the  simplest  case  for  a  class  of  equations  that  have 

been  called  "A-w  -systems"  [10],  [11],  [15],  [22],  [23].  We  now  try  to  solve 

LV 

(6.5)  by  the  ansatz  *  re  ,  r  *  C  sech  kx,  *  B  tanh  kx.  After  some 
elementary  calculations  this  leads  to  the  equations 

_  _  A 

k  -  B  *  ±  Re  a1  +  u  Re  a3 

(6.6)  -  2k2  +  B2  *  C2Re  a2 

A 

-  2Bk  »  i  In  a^  +  u  Im  a^ 

2  '• 

3Bk  *  C  Im  a2 

From  the  fourth  and  second  equation  of  (6.6)  we  find 

3Bk  Re  a2  +  (2k2  -  B2)  In  #2  »  0 

A 

which  can  be  solved  by  B  *  Ak,  where  (provided  Ira  a2  *  0) 
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According  to  the  fourth  equation  of  (6.6),  X  must  have  the  same  sign  as 

■V 

1m  a2 ,  which  is  achieved  by  choosing  the  minus  sign  in  the  numerator . 

We  now  insert  B  *  Xk  into  the  first  and  third  equation  of  (6.6),  thus 
obtaining 

k2(1  -  X2)  =*  ±  Re  a^  +  u  Re 

-2Xk2  =  ±  Im  a^  +  w  1m  a3 

2  2 

»»>  k  {(1  -  X  )  Im  a3  +  2X  Re  a3)  =  ±(Re  a^  Im  -  Im  a^  Re  a3) 

The  right  side  of  this  last  equation  is  not  zero,  according  to  what  we  have 
said  about  a^  and  a^  above.  We  can  resolve  the  equations  with  respect 

a 

to  k  and  o>  if  the  following  holds 

(vi)  (1  -  X2)  im  a3  +  2X  Re  a3  *  0 

If  ft  »  0,  (6.6)  can  be  solved  by  B  *  0,  provided  that  the  following 

holds: 

(vii)  If  ft  *  0,  then  Re  is  negative  and  Im  a3  ^  0 
We  find  thus 


u(x,t)  '  uq(m)  +  2C(Ce 


+  c.c.) 


+k£x  i(t-BEx) 
e 

This  means  that  asymptotically  we  have  exponentially  damped  waves, 
propogating  in  opposite  directions  for  x  +  ±®. 

ii p 

A  second  ansatz  for  a  solution  of  (6.6)  is  ■  re  ,  r  *  C  tanh  kx, 

V'  ”  B  tanh  kx.  This  ansatz  leads  to 

2 

-2k  *  ±Re  +  w  Re  a^ 

A 

(6.7)  3kB  »  i  In  a(  +  u  In  a^ 

2  2  2  ^ 

2k  -  B  *  C  Se  a2 

2  ' 

-3kB  ■  C  In  >2 

The  last  two  equations  can  again  be  solved  by  the  ansatz  B  *  Ak  with  the 

same  expression  for  A  as  before,  this  time,  however,  the  plus  sign  must  be 

chosen  in  front  of  the  square  root.  The  first  two  equations  now  lead  to 

2 

-2k  =■  t  Re  a,  +  hi  Re  a, 

1  3 

2  * 

3k  A  -  ±  Ira  a1  +  <■»  Im  a3 

2 

**>  k  (-2  IB  »3  -  3 ARe  a3>  *  ±  (Re  a^Im  a^  -  Im  ai  Re  a3) 

These  equations  can  be  resolved  w.r.  to  k  and  u  if 

(vi) ’  21m  a3  +  3A  Re  a3  *  0  . 

For  Im  a2  =  0,  once  can  again  resolve  (6.7)  by  B  =  0,  provided  that 

(vii) '  If  In  *  0,  then  Re  a2  is  positive  (note  that  this  is  the 

opposite  of  condition  (vii)),  and  1m  a3  ?  0. 

We  thus  find 
Proposition  6.3: 

If  (vi)’  or  (vii)'  respectively  are  satisfied,  then  the  ansatz 

i  <fi  * 

-  re  ,  r  *  C  tanh  kx,  '  ■  B  tanh  kx  yields  a  solution  (x)  to  (6.5) 

This  solution  is  an  odd  function  of  x  and  converges  to  periodic 

wave  trains  propagating  in  opposite  directions  as  x  ♦  ±*».  It  agrees  with  the 

solution  considered  in  [38]  as  a  model  for  one- dimensional  spiral  waves.  In 
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that  paper,  (6.5)  and  its  analogue  in  two  space  dimensions  are  considered  for 

* 

the  special  case  Re  =  Ira  a1  *  0.  The  solution  (x)  is  also  the  simplest 
of  the  solutions  the  existence  of  which  has  been  proved  in  [23]  for  the  case 

■v 

that  Im  a2  is  small. 

we  now  want  to  prove  that,  for  sufficiently  small  e  ,  there  exists  a 

0  * 

solution  having  the  same  characteristics  as  and  a  ,  respectively,  we 

0  A 
begin  with  the  case  of  .  Linearizing  (6.4)  with  respect  to  o  and  o>,  we 

obtain  the  following  inhomogeneous  problem  that  has  to  be  discussed 

b'  -  0  =  f 

1  2  1 


(6.8) 


S2  *  aiB1  *  a2(°?>2l1  “  2a2  °X8r  *3“* 1  '  a 3 "  f 2 


for  which  we  write  briefly. 


L0  +  (0,-a3aJ)S)  *  f 


As  in  chapter  3,  we  want  to  study  properties  of  this  linearized  operator  as  a 


mapping  from 


R  x  Vi  =  {(fW  e  cb+1(R'c2)l  Br>{x)'  ife.  BIk)(x) 

exist  for  0<k<m+1 ,  0^  (x)  =*  0^ (-x) ,  B^tx)  *  -  02<-x)} 


Fm  =  { (f  1  ,f 2 )  e  C™(R,C2)|1^j  fjk),  ^m<»  f{k)(x)  exist  for  0<k<m, 
f  ^  (x)  =  -f^-x),  f2(x)  *  f2(-x)} 

2 

The  point  0  e  C  is  a  saddle  point  for  (6.4),  which  has  two  stable  and  two 

unstable  directions .  Hence  the  same  arguments  as  in  §  3  show  that  for  any 

given  f  e  F  and  Si  e  R,  we  can  find  a  bounded  solution  0  on  [0,“),  and 
m 

two  initial  conditions  at  0  are  left  arbitrary.  These  initial  conditions 


can  be  matched  by  adding  multiples  of  the  solutions  of  the  homogeneous 
problem,  which  are  given  by  (c^  ,a  )  and  (i<*^,  i°2)  •  As  in  §  3,  the 
solution  on  [0,°»)  extends  to  a  solution  on  all  of  R,  if  02 ( 0 )  is  zero. 


O'  0 

As  we  can  verify  from  (6.6),  Re  (0)  *  0,  but  i<»2(x)  vanishes  at  zero. 

This  means  that  by  appropriate  choice  of  the  elements  of  the  nullspace  we  can 

adjust  one  initial  condition  at  x  ■  0,  the  other  must  be  matched  by 

appropriate  choice  of  fl.  In  other  words,  L  (as  an  operator  from  U 

ni+i 

into  F  )  has  a  one  dimensional  nullspace  and  the  range  has  codimension  one. 
m 

We  assume 

(viii)  (0,-a3a°)  is  not  in  the  range  of  L 
Now  the  same  arguments  as  in  $  4  lead  to  the  result 
Theorem  6.4? 

Assume,  conditions  (i)-(viii)  hold.  Then  for  each  e  in  a  neighborhood 

A 

of  0  there  exists  o>  e  R  for  which  (6.4)  has  a  one-parameter  family  of  non¬ 
vanishing  solutions,  which  are  even  in  x  and  approach  0  as  x  ±°°.  The 
solutions  in  this  one  parameter  family  differ  from  each  other  only  by  a  shift 
in  the  time  variable . 

Remark ; 

Condition  (viii)  is  of  a  "generic"  type  in  the  sense  that  it  requires  a 
certain  quantity  not  to  vanish .  We  have  not  succeeded  in  giving  an  explicit 
criterion,  when  (viii)  is  true.  In  one  particular  case,  however,  this  can 
easily  be  seen,  namely,  assume  that  Im  a^  »  0,  and  Im  a^  *  0.  Then  it  can 
be  shown  that  the  functional 

00 

(f1,f2)  ♦  im  /  ( — (x)f1 (x)  +  a°(x)f2(x))dx 

«QO 

annihilates  the  range  of  L,  and  clearly,  (O^a^)  is  not  in  the  nullspace 
of  this  functional .  It  would  be  interesting  to  know  whether  there  are 

parameter  values  for  which  (viii)  is  false. 

* 

For  the  case  of  the  solution  c^,  we  have  to  introduce  some  new 
definitions  of  spaces,  since  we  are  now  dealing  with  solutions  approaching  a 
periodic  limit  at  infinity  rather  than  a  constant,  and,  moreover,  the  symmetry 
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ks 


properties  are  different  (a^  is  odd  rather  than  even),  which  means  we  shall 
have  to  make  use  of  the  reversibility  under  R  rather  than  under  R. 

■v 

Throughout  the  following  discussion,  we  shall  assume  that  Im  *  0,  whence 
B  *  0  ,  and  it  is  no  restriction  to  assume  B  >  0,  (otherwise  change  the 
signs  of  both  B  and  k,  which  just  corresponds  to  a  symmetry 
transformation) . 

Definition  6«5: 

Let  o  be  a  positive  real  number.  Then  Z?  (Rn)  denotes  the  space  of 

*  ,o 

all  functions  (u(x),  v(x))  :  R  ♦  ^  (Rn)  x  A1  «  (Rn)  such  that  the  following 

m  m-  /2 


1 .  u  and  v  respectively  are  C  -functions  of  x  e  R 

2.  On  (0,®)  u  and  v  respectively  can  be  represented  as  the  sum  of  a  2it- 

£  —Ox 

periodic  C  -function  and  the  product  of  e  with  a  function,  whose 

first  £  derivatives  are  continuous  and  bounded  and  converge  to  zero  at 

infinity  . 

3 .  All  odd  Fourier  components  of  u  (with  respect  to  t)  are  odd  functions 
of  x,  all  the  even  Fourier  components  are  even  functions  of  x,  and  vice 
versa  for  v . 

a 

Z.  is  defined  in  the  same  way,  but  condition  3  is  to  be  taken  the  other  way 

J6-  r  U 

round.  Clearly  there  is  a  natural  choice  for  a  norm  in  z?  ,  and  we  omit  a 

l, a 

detailed  definition. 

Me  fix  a  <J  >  0,  which  is  to  be  chosen  sufficiently  snail,  and  we  look 

for  solutions  to  (6.4)  in  Z^™^  which  lie  in  a  neighborhood  of 
*  *  ^ 

(aj(x),  <*2  (x) ,  z  *  0)  .  Since  we  wish  to  fix  the  period  of  the  periodic  part 
to  2*,  a  scaling  factor  Y  must  be  introduced,  for  e  -  0  this  factor 
equals  B  (i.e.  in  (6.4)  we  put  x  *  Yx  ,  but  we  shall  again  write  x  for 

>■  *  h 

x)  .  Again  we  must  discuss  the  linearization  of  (6.4)  at  (a,^ ,  a^,  0).  The 
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third  equation  causes  no  problems,  as  can  be  seen  from  the  results  of  §  4  and 
the  fact  that 


<6Y  -  A ( 0 ) )  (e”0Xf (x) )  -  e“aX(ey  -  eyo  -  A(0))f(x) 

ax  ax 

Let  us  now  investigate  the  linearization  of  the  first  two  equations  of 

*  * 

(6.4)  at  (al  (x)  #  a2(x)|,  i  *e  •  we  have  to  discuss  the  nullspace  and 

2  in  A 

oodimension  of  the  operator  D:D(D)  C  r  x  Z?  +  2^  given  by 


★  i 


*  i 


(6.9)  D(fl.r#61#62)  =»  (Ta1  +  -  02,  r<»2  +  b$2  *  a  ^ 


-■3“6,  -  2*2  V>,  -  *2  "V2®1  '  *J°V 


Let  f  -  (f  ,f  )  e  2-  be  given  and  consider  the  problem 

»  2  f  0 

=  f.  Due  to  symmetry  properties  it  is  again  enough  to  find 

solutions  8  on  the  half-line  [0,“),  which  satisfy  8^(0)  =0.  we  split 

* 

8,f  and  a  into  the  periodic  and  exponentially  decaying  parts,  denoting 
them  by  indices  p  and  d.  The  periodic  part  of  (6.9)  reads 


•  * '  - 
b84  -  B„  -  f,  -  To  =:  f. 

Ip  2p  Ip  Ip  Ip 

(6.,0)  BS2-  *  .,»)p  -  .,*1p  -  -  2.'2l»;pl261p 

*1  ft  % 

*  f_  -  To  +  a,8a  = :  f 
2p  2p  3  Ip  2p 


*  +ix  i# 

Here  we  have  a  *  ±Ce  e  ,  where  C  is  the  constant  named  so  in 
>P 

proposition  6.3,  the  +  or  -  sign  agrees  with  the  sign  of  k,  and 
00 

*  =  B  /  (tanh  (kx)  -  1)dx  for  k  >  0 
0 


or 


*  *  B  /  (tanh  kx  +  1 )  dx 
0 


for  k  <  0 


Proposition  6.6: 


The  linear  operator  represented  by  the  left  side  of  (6.10)  (in  the  space 


of  2n-periodic  C  -functions)  has  a  one  dimensional  nullspace  spanned  by 

•*  *  *  • 

(a  ,  a  )  and  oodimension  one.  The  functional  annihilating  the  range  is 

Ip  2p 

given  by 


tm 

(frf2)  +  im  / 


2*  ?2ix  -2i® 

- - - -  (a  f  +  at  ) 

0  a  2P  1  ip  i 


Proof  s 


+ix 

In  (6.10)  we  put  3  -  e  e  3,  thus  obtaining 


(6.11) 


\  I  >•  V  X<  V  «4  4% 

bB  ±  IbB,  -  B„  -  e  e  1 V 

Ip  Ip  2p  Ip 


b3  '  ±  iB3_  +  al.  -  awB.  -  a.  |ot  |2I  -  2a,  |o.  |23. 

2p  2p  1  Ip  3  Ip  2  Ip  Ip  2  Ip  Ip 

;  +ix  -i® 

*  f-  e  e 

2p 


In  order  to  find  the  nullspace,  we  make  the  ansatz 

i  inx  -inx 

3.  =  Pe  +  oe 

Ip 

This  leads  to  the  equations 

2  _  2  "2  '  2— 

( (— n  +  2n)B  -  a2C  )P  -  a2C^J  *  0 

2  2  %  2  »  2— 

((-n  ±  2n)B  -  a2C4)0  -  a£Tp  =  0 

Non-trivial  solutions  exist  only  if 

o  O  >-  2  2  2^2  %  >  A 

((-n  ±  2n)B  -  a2C  )((-n  *  2n)B  -  a2C  )  *  a2a2C 

4  24  222  2  2  ' 

<==>  (n  -  4n  )B  +  2n  B  C  Re2a  ±  4iB  C  nlm2a  *  0 

==>  n  =  0 

From  this  it  is  easily  seen  that  the  nullspace  is  one  dimensional.  Since  the 
resolvent  is  compact,  so  is  the  oodimension.  One  can  easily  verify  explicitly 


that  the  functional  given  above  annihilates  the  range 


I 


We  now  turn  to  the  discussion  of  the  decaying  part.  It  is  determined  by 
the  equations 

•  *i  * 

BB1d  “  02d  “  f1d  ”  r°1d  "!  f Id 
(6.12)  B«Jd  *  -  ^V^ld  -  2^8,,, 

"  f2d  -  T*U  +  a3fl<,*d  *  *2  <Vftp  *  “j'hVd  S1p  -  f2d 

Proposition  6.7; 

A  A 

For  any  given  right  side  (f,  .,  fn.)t  equation  (6.12)  has  a  one- 

id  2d 

parameter  family  of  solutions  on  [0,“),  which  approach  zero  exponentially  as 
x  ♦  •. 

Proof ; 

it 

In  the  limit  x  ♦  “,  all  terms  on  the  left  of  (6.12)  that  contain  a„_. 

la 

can  be  regarded  as  a  perturbation  (cf.  §  3).  Hence  it  suffices  to  show  the 

jbix 

proposition  is  true  if  we  drop  these  terms.  Again  we  substitute  3  -  e  e  3, 

and  we  are  left  with  the  same  left  hand  side  as  in  (6.11),  except  that  now  we 

— Ox 

must  look  for  solution  decaying  to  zero  at  a  rate  of  e  rather  than  for 
periodic  solutions. 

Going  through  the  same  steps  as  we  did  following  (6.11),  we  obtain  the 

following  equation  for  the  characteristic  exponents  belonging  to  (6.10): 

44  42  222^  22  % 

A  b  +  4B  A  -  2A  B  C  Re  a^  ♦  4B  C  A  Im  a^  *  0 

One  eigenvalue  is  A  =■  0,  and  it  is  simple  unless  Im  a^  *  0,  which  we  have 

excluded.  Using  (6.7),  we  obtain  for  the  remaining  eigenvalues: 

AV  +  6Ab2  -  4Ak2  ±  12kB  -  0 

which  is  solved  by  A  *  +  2—,  leaving  the  following  equation  for  the  remaining 

B 


eigenvalues 


A2b2  *  2Akfl  +  6B2 


t2kB  ±  l/4k2B2  -  24B4 


The  eigenvalue  T  2—  is  negative,  and  the  last  two  eigenvalues  have  positive 
B 
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real  parts.  This  yields  the  proposition. 

Solutions  on  [0,*)  yield  solutions  to  (6.12)  on  all  of  R  if 

^(0)  =  0.  On  (O,00)  we  have  two  linearly  independent  solutions  of  the 

*  *  *  •  *  •  *  • 

homogeneous  problem,  given  by  (ia^ia^)  and  (a^  ,ct  •  Since  (0)  *  0, 
but  a*(0)  *  0,  addition  of  a  solution  of  the  homogeneous  problem  cam  only  be 
used  to  match  one  of  the  two  initial  conditions.  Hence  we  see: 

■x. 

The  operator  D  =  has  a  one-dimensional  nullspace  and  its  range 

has  oodimension  2  (one  coming  from  the  condition  8  (d)  =  0  and  one  from  the 

codimension  for  the  periodic  part)  .  Thus  D  is  onto  iff 

*  •  * «  *  • 

(viii)'  The  vectors  (<*  ,a  )  and  (0,-a  a  )  span  a  complement  to  the  range 

1  £  <3  1 

of  D. 

We  do  not  know  how  to  check  (viii)'  explicitly,  tut  one  can  expect  it  to  hold 
for  almost  all  parameter  values. 

Theorem  6.8: 

Suppose  (i)-(v) , (vi) ' , (viii) '  hold  and  Im  *  0.  Then  in  a 
neighborhood  of  e  =  0  there  exist  Y(£)»  w(e),  for  which  (6.4)  (with  x 
scaled  by  Y)  has  a  one-parameter  family  of  non-periodic  solution  that  lie  in 
0*  The  solutions  in  this  one-parameter  family  are  again  distinguished 
only  by  a  time  shift. 
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7 .  Stability 

This  chapter  deals  with  the  question  whether  the  solution  provided  by 
theorem  6.4  is  stable.  We  shall  prove 
Theorem  7 . 1 : 

If  v  is  real  (e.g.,  if  all  diffusion  coefficients  are  equal),  Im  a2 

■s* 

is  small  and  Re  a^  <  0,  the  solution  given  by  theorem  6.4  is  unstable. 

Proof : 

It  is  convenient  to  introduce  artificially  a  second  time  variable  t. 

I.  e.,  we  proceed  as  follows:  We  seek  solutions  to  (6.1),  which  have  the  form 

3 

u(t  +  t)  and  are  periodic  in  t.  The  operator  must  then  be  replaced  by 

3  3 

.  If  we  then  go  through  the  same  transformations  that  led  to  (6.4), 
we  arrive  at  the  system 


(7.11  <4  -  *  V)  *  (b*D,D-V|f 

+  0( |e|  +  »£») 

*  _  3a 

7-  *  -v  n  v*b  D  -*i  3z  1  1 

ez  -  a(o)z  +  (i  -  ■  -  )d  ’(-g|  +  7  aT1  *  y)  +  o(e) 

(b  D,y) 


Let  us  now  assume  that  v  is  real  and  a2  is  real  (and  negative) .  As  a 

A 

consequence,  a^  is  imaginary,  a^  +  is  real  and  the  plus  sign  must  be 

chosen  in  front  of  a1  (cf.  (6.6))  for  obtaining  the  solution  a®(x).  we 

0  0 

linearize  (7.1)  at  the  solution  ai  -  a^(x),  a2  ■  ®2(x),  z  *  0,  and  seek 

c2Xt 

solutions  of  the  linear  equation  which  are  proportional  to  e  .  For  e  *  0 
this  yields  the  following  eigenvalue  problem  for  X 

-  *1S1  +  “2VaV  +  Vl'-l'2  ‘  V81  *7®, 

When  31  is  restricted  to  be  real  valued,  we  obtain 


(7.2) 


+  V  +  3*2<al)2)01 
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0 '  2 
The  derivative  is  a  solution  for  X  -  0.  The  operator  in  h 

represented  by  the  right  hand  side  of  (7.2)  is  self-adjoint,  and  its 

eigenfunctions  are  critical  points  of  the  functional 

CO 

e  9  2  *  *022 

J  -  P1  -  (a^  a30J  +  3a2(a“r)e‘dx 

-00 

2 

on  the  unit  sphere  of  L  .  It  is  easy  to  prove  that  the  maximizing  function 

does  not  change  its  sign  (this  is  a  well-known  principle  in  quantum  mechanics) 

O' 

and  thus  cannot  be  a  .  Hence  there  is  an  eigenvalue  X  >  0,  and  the 
linearized  equation  has  an  exponentially  growing  solution .  Standard 
perturbation  theory  shows  that  this  property  is  preserved  under  small 

■v 

perturbations  in  e  and  a2,  which  gives  the  theorem. 
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